Edith Cowan University 
Research Online 


Research outputs 2022 to 2026 


4-8-2022 


Assessment of wettability and rock-fluid interfacial tension of 
caprock: Implications for hydrogen and carbon dioxide geo- 
storage 


Muhammad Ali 
Bin Pan 

Nurudeen Yekeen 
Sarmad Al-Anssari 


Amer Al-Anazi 


See next page for additional authors 


Follow this and additional works at: https://ro.ecu.edu.au/ecuworks2022-2026 


G Part of the Chemical Engineering Commons 


10.1016/j.ijhydene.2022.02.149 

Ali, M., Pan, B., Yekeen, N., Al-Anssari, S., Al-Anazi, A., Keshavarz, A., ... & Hoteit, H. (2022). Assessment of 
wettability and rock-fluid interfacial tension of caprock: Implications for hydrogen and carbon dioxide geo-storage. 
International Journal of Hydrogen Energy, 47(30), 14104-14120. https://doi.org/10.1016/j.ijhydene.2022.02.149 
This Journal Article is posted at Research Online. 

https://ro.ecu.edu.au/ecuworks2022-2026/534 


Authors 
Muhammad Ali, Bin Pan, Nurudeen Yekeen, Sarmad Al-Anssari, Amer Al-Anazi, Alireza Keshavarz, Stefan 
Iglauer, and Hussein Hoteit 


This journal article is available at Research Online: https://ro.ecu.edu.au/ecuworks2022-2026/534 


INTERNATIONAL JOURNAL OF HYDROGEN ENERGY 47 (2022) 14104-14120 


Available online at www.sciencedirect.com International Journaljof 
HYDROGEN 


ScienceDirect ENERGY _ 


ee at? = 
ELSEVIER journal homepage: www.elsevier.com/locate/he 


S 


Assessment of wettability and rock-fluid interfacial 
tension of caprock: Implications for hydrogen and 
carbon dioxide geo-storage 


Check for 
updates 


Muhammad Ali **, Bin Pan’, Nurudeen Yekeen ©, Sarmad Al-Anssari “, 
Amer Al-Anazi °, Alireza Keshavarz °, Stefan Iglauer °, Hussein Hoteit ” 


* Physical Science and Engineering Division, King Abdullah University of Science and Technology (KAUST), Thuwal, 
23955, Saudi Arabia 

P School of Civil and Resource Engineering, University of Science and Technology Beijing, Xueyuan Road, Beijing, China 
© Department of Chemical & Petroleum Engineering, Faculty of Engineering, Technology and Built Environment, UCSI 
University, 56000, Kuala Lumpur, Malaysia 

4 Department of Chemical Engineering, University of Baghdad, Baghdad, 10071, Iraq 

© School of Engineering, Edith Cowan University, Joondalup, WA, 6027, Australia 


HIGHLIGHTS GRAPHICAL ABSTRACT 


e The equilibrium contact angle (g) 
increases with pressure and 
organic concentration but de- 
creases with temperature. 

e The solid-gas interfacial tension 
(Ysg) decreases with pressure and 


Caprock 


organic concentration but in- 
Hydrogen 


creases with temperature. on, Ges = 
e The solid-liquid interfacial tension 

(ys) and 0g decreases, whereas ysg 

increases with increasing alumina 

nanofluid concentration. 
e Mica hydrogen wettability is less 

than mica CO, wettability in all 


physio-thermal conditions. 
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(6g) and solid-liquid and solid-gas interfacial tensions (ys, and ysg). However, due to the 
technical constraints of experimentally determining these parameters, they are often calculated 
based on advancingand receding contact angle values. There is a scarcity of 6g, ys, and ysg data, 
particularly related to the hydrogen structural sealing potential of caprock, which is unavailable 
in the literature. Young's equation and Neumann's equation of state were combined in this study 
to theoretically compute these three parameters (0g, Ys, and ysg) at reservoir conditions for the 
Hz and CO, geo-storage potential. Pure mica, organic-aged mica, and alumina nano-aged mica 
substrates were investigated to explore the conditions for rock wetting phenomena and the 
sealing potential of caprock. The results reveal that 6g increases while ysg decreases with 
increasing pressure, organic acid concentration, and alkyl chain length. However, ysg decreases 
with increasing temperatures for H3 gas, and vice versa for CO2. In addition, #¢ and ys, decrease, 
whereas ysg increases with increasing alumina nanofluid concentration from 0.05 to 0.25 wt%. 
Conversely, Og and ys, increase, whereas ysg decreases with increasing alumina nanofluid 
concentration from 0.25 to 0.75 wt%. The hydrogen wettability of mica (a proxy of caprock) was 
generally less than the CO, wettability of mica at similar physio-thermal conditions. The 
interfacial data reported in this study are crucial for predicting caprock wettability alterations 
and the resulting structural sealing capacity for UHS and CGS. 
© 2022 The Author(s). Published by Elsevier Ltd on behalf of Hydrogen Energy Publications 

LLC. This is an open access article under the CC BY license (http://creativecommons.org/ 


licenses/by/4.0/). 


Introduction 


The rapid industrialization and increasing world population 
(at a geometrical progression rate) have increased global en- 
ergy demand and consumption [1]. Fossil fuel has remained 
the principal global energy source for several decades [2,3], but 
the world reserve is depleting faster than the discovery rate 
[4,5]. Moreover, whenever fossil fuels are burned, they release 
large amounts of greenhouse gases into the atmosphere, 
contributing to global warming [6—10]. 

Underground hydrogen (H2) storage (UHS) and carbon di- 
oxide (CO2) geo-storage (CGS) have been identified as two 
effective methods of providing sufficient global energy needs 
and creating a low-carbon global economy by reducing carbon 
emissions [11—15]. The primary aims of UHS and CGS are 
substituting fossil fuel with decarbonized and clean fuels (H3) 
[14,16,17], converting excess renewable energy into H through 
electrolysis [14,18], and capturing CO, from stationary sources 
or the atmosphere and sequestering it in geologic formations 
[19—21]. In addition, CO, sequestration involves the permanent 
and secure storage of the captured gas in underground for- 
mations. In contrast, UHS involves the storage of H3 in geo- 
storage locations (such as depleted oil and gas reservoirs, 
coal seams, deep saline aquifers, and organic-rich shales) and 
its withdrawal at the period of peak demand [14,16,17,22—24]. 

However, the success of UHS and CGS is greatly influenced 
by several thermo-physical and petro-physical factors [18]: 
the wettability of the rock-brine-gas systems and the inter- 
facial tension (IFT) between the rock and fluids (i.e., gases or 
brine) [25,26]. These crucial parameters govern the fluid dis- 
tribution in the reservoir, the injection of CO; and H; into the 
reservoir, and the subsequent H, withdrawal [27—32]. Evalu- 
ation of the rock wetting phenomenon and rock-fluid IFT is 
crucial to assess the gas storage capabilities of the storage- 


rock and caprock formations. It also provides operational 
constraints to avoid leakage hazards and ensure containment 
safety at realistic storage conditions [27,33—35]. 

Generally, reasonable efforts have been invested in 
measuring the contact angle of rock-brine-CO, systems in 
previous studies [25,36—38]. However, only four studies have 
reported contact angle datasets for rock-brine-H2 systems 
[26,34,39,40]. Measuring the contact angle of rock-brine-H, 
systems is challenging due to high H, reactivity, volatility, 
compressibility, and metal embrittlement [34,41]. 

The possibility of industrial-scale Hz storage in sandstone 
reservoirs (mainly consisting of quartz) with the presence of 
organic acid contamination was investigated by Iglauer et al. 
[31] and Ali et al. [34]. Iglauer et al. [31] demonstrated that 
higher contact angles were measured at an elevated temper- 
ature and pressure with the presence of a minute concentra- 
tion of organic acid. The results from Ali et al. [34] agreed with 
those from Iglauer et al. [31], confirming that the intrinsically 
strongly water-wet clean quartz substrate became H2-wet (the 
advancing contact angle increased from 40.8° to 91.3°) when 
aged in lignoceric acid (107° M) at 25 MPa and 323 K. Similarly, 
Ali et al. [26] found that pure mica substrates became 
intermediate-wet at higher pressures and lower tempera- 
tures, whereas organic acids converted them into H,-wet. 

In contrast, forrock-fluid IFT, only gas-brine IFT (Ygas-brine) can 
be conveniently measured in the laboratory [27]. Rock-gas IFT 
data (Yrock_-H, aNd Yrock-co,) are scarce in the literature, as these 
parameters cannot be determined by experimental estimations 
[27,42,43]. For rock-brine-CO, systems, researchers computed 
these data theoretically through the combination of Neumann's 
equation of state and Young's Laplace equation [27,29,30,42—47]. 
Arif et al. [43] found that mica-CO, IFT (Ymica-co,) Was higher at 
high temperatures and lower at high pressures. Similarly, Pan 
etal. [46] demonstrated that yquarz-co, is considerably lower than 
Yquariz-cH, at similar thermo-physical conditions. 
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To the best of our knowledge, only Pan et al. [27] computed 
rock-fluid IFT for rock-brine-H, systems using a semi- 
empirical technique, suggesting that evaluating rock-fluid 
IFT data for rock-H>-brine systems at geo-storage conditions 
requires further investigation. Pan et al. [27] reported that the 
values of rock-gas IFT were smaller at a higher temperature 
and elevated pressure and were reduced with an increased 
organic acid concentration. However, the reported rock-fluid 
IFT data were primarily for quartz, a representative mineral 
in a sandstone formation. Thus far, no data are available for 
rock-fluid IFT for caprock representative minerals (i.e., mica). 
Hence, the associated structural trapping capacities of the 
caprock during Hz and CO, storage conditions are still un- 
known, which motivated this work. 

Initially, the advancing (6,) and receding (6,) contact 
angles of the mica-H2-brine and mica—CO,—brine systems 
were measured for pure, organic-aged, and nano-aged sub- 
strates. Afterward, we computed the mica-fluid IFT data 
(Ymica—Ho, Ymica—COs, and Ymica—brine ) Using the combination of 
Young's equation and Neumann's equation of state in geo- 
storage conditions. 

The focus has been on investigating the influence of tem- 
perature, pressure, gas type (CO, and H3), organic acid con- 
centration, alkyl chain length, and nanofluid concentration on 
contact angles and rock-fluid IFT data in mica-H2-brine and 
mica—CO,—brine systems. To our knowledge, this study pro- 
vides the first reported rock-fluid IFT data for H, and a caprock 
representative mineral (mica) at subsurface conditions. The 
data indicate favorable conditions for effective structural 
trapping capacities of the caprock and gas containment se- 
curity during UHS and CGS. 


Materials and methods 
Materials 


The details of the materials used for this study are presented 
in Tables 1 and 2, where mica (muscovite) is used as a 


Table 1 — Materials used for this research. 


Materials Suppliers Other details 


Pristine mica Ward's natural LxWxH=20x15x 3mm 


science 
NaCl Chemlab >99.9 mol % purity 
Deionized water David Gray 0.02 mS/cm electrical 


conductivity 


N-decane Chem-supply 99.9 mol % pure 
Nitrogen BOC 99.999 wt% ultra-pure (gas 
code 234) 
CO, BOC 99.999 wt% ultra-pure (gas 
code 082) 
H2 BOC 99.999 wt% ultra-pure (gas 
code 240) 
Hydrochloric Sigma Aldrich >99% pure, 37.5% 
acid concentration 
Hydrophilic Sigma Aldrich Particle size (13 nm), surface 
aluminum area (85—115 m?/g), purity 


oxide (Al,03) 
nanoparticles 


(99.8 wt%) 


representative mineral for the caprock. The contact angle 
measurements for H, and CO, were conducted with a 10 wt% 
NaCl aqueous solution at reservoir conditions. Four types of 
organic acids with varying carbon chain lengths (lignoceric 
acid C24, stearic acid Cy, lauric acid C13, and hexanoic acid Ce) 
were used as representative organic molecules in geological 
formations to age the mica substrates [35,48—51] (Table 2). The 
mica samples were cleaned with ultra-pure nitrogen before 
aging to eliminate impurities from the mica surface. The 
various organic acid concentrations were prepared using n- 
decane as the base solution. The mica substrate surfaces were 
cleaned with methanol, acetone, and toluene (99.9 mol % 
pure, Rowe Scientific). Hydrophilic aluminum oxide (Al)03) 
nanoparticles were used to improve the wetting phenomenon 
of organic-aged mica substrates. 


Cleaning and aging procedure for mica substrates with 
organic acids 


Pure mica substrates must be cleaned because they contain 
impurities in the initial condition that could influence the 
experimental results. Therefore, mica substrates were first 
cleaned with deionized water, and then the samples were 
blown with ultra-clean nitrogen and dried at 353 K for 30 min 
to eliminate the thin water film. Afterward, the residual 
minute organic impurities were removed by treating the mica 
samples with air plasma (Diemer Yocto instrument) for 20 min 
[35,52,53]. 

Realistic caprock conditions were replicated through the 
following process. First, mica surfaces were inserted into 
brine (2 wt% NaCl) for half an hour at room conditions to 
ionize the surface. A pH of 4 in the brine solution was main- 
tained by adding HCl droplets. This procedure promotes the 
adsorption of the organic acid on the mica surfaces to mimic 
realistic geological conditions [35,49,51,54]. 

Subsequently, ultra-pure nitrogen flow was used to remove 
the residual brine from the modified mica surfaces. Then, the 
HCl-brine ionized substrates were put in diverse solutions of 
n-decane or organic acid mixtures, with concentrations of 
10-?, 10-3, 10-°, 1077, and 107°? mol/L. The mica surfaces were 
aged in organic solutions for 1 week to imitate exposure to 
formation water for prolonged geological stints [35,49,51,54]. 
Next, mica substrates were cleaned with ultra-pure nitrogen 
followed by drying for 24 h at 343 K to remove any excessive 
organic or n-decane solution before contact angle measure- 
ments. The mica substrates became hydrophobic after aging 
in the organic solutions because of the chemisorption of 
organic acids on the mica surfaces [35,55]. 


Mica substrate characterization 


Measurement of the total organic carbon 

The total organic carbon (TOC) of the aged and unaged mica 
substrates was measured using the pyrolysis method with the 
Rock-Eval 6 instrument (Vinci Technologies) [35]. The TOC of 
the clean mica substrates was insignificant (around 500 ppm 
due to air contamination). However, when the mica substrate 
was aged in organic acids, the TOC increased up to 1800 ppm 
owing to organic acid chemisorption. The TOC results 
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Table 2 — Organic acid details. 


Name 

Lignoceric acid Cy4H4g0 368.63 
Stearic acid CygH3602 284.4772 
Lauric acid Cy2H 2402 200.318 
Hexanoic acid C6H1202 116.158 


Formula Molecular mass (g/mol) Number of carbon atoms pH (pKa) 


Physical state Chemical structure 


24 7.4 Solid (0) 


a Sek 


18 6 Solid (0) 


a H 


12 BES) Solid (6) 


Ree ae 


4 Liquid (0) 
ienen a 


presented in Fig. 1 are responsible for altering the wetting 
characteristics of the mica substrates. 


Surface characterization of mica substrates 

Energy dispersive spectroscopy coupled with a field-emission 
scanning electron microscope (SEM; Oxford Instruments) and 
atomic force microscopy (Nanosurf, Flex-Axiom, and 
Controller C3000) were used to analyze the elemental 
composition and surface topography of the clean and aged 
mica substrates. The pure mica surfaces were very smooth, 
and their surface roughness was 1—2 nm. Moreover, organic 
acid adsorption has very little effect on the surface roughness 
(and, thus, the contact angle), and the root mean square 
ranges from 200 to 460 nm. Many wettability studies have 
reported that a surface roughness of less than 1 um does not 
significantly influence the contact angle measurements 
[56—58]. 

The SEM micrographs of the mica substrate before and 
after aging in stearic acid (107° mol/L concentration) are pre- 
sented in Fig. 2, exhibiting a uniform coverage of the organic 
acid on the mica surface. A significant increase in the carbon 
concentration on the mica surface was observed for the aged 
mica substrates. For instance, the mean carbon content of the 
pure mica increased from 3.2 + 0.5 to 5.0 + 1.6 wt% when aged 
in lignoceric acid solution, from 3.2 + 0.6 to 4.8 + 1.6 wt% when 
aged in stearic acid solution, from 3.1 + 0.7 to 4.4 + 1.7 wt% 
when aged in lauric acid solution, and from 2.9 + 0.7 to 
3.9 + 0.8 wt% when aged in hexanoic acid solution (Table 3, the 
total energy dispersive spectroscopy values of each mica 
organic added substrate is depicted in the supplementary file 
in Table S1). These results can be attributed to the organic acid 
chemisorption on mica surfaces [34]. 


Preparation of the alumina nanofluid 


Different nano-alumina concentrations (0.05—0.75 wt%) were 
prepared using the sonication technique to explore their 


capacity to alter the wetting phenomenon of the mica—COp. 
—brine system. The desired Al,O3 nanoparticle concentrations 
were dispersed in ionized water using a magnetic stirrer to 
ensure homogeneity. Afterward, the solutions were sonicated 
using an ultrasonic homogenizer (Sonics and Materials 
Incorporation, USA) for 15 min at a frequency of 20 kHz with 
40% sonication amplitude and an energy of 9500 J [59]. 


Determination of nanofluid stability 

The stability of the prepared nanofluids dispersions must be 
ascertained because agglomerations of nanoparticles in 
deionized water or brines (associated with nanoparticle hy- 
drophobicity and surface energy) could decrease the wetta- 
bility reversal properties [50,60,61]. The stability and phase 
behaviors were estimated from the hydrodynamic particle 
sizes and zeta-potential values measured using a Malvern 
nano zeta-sizer (Model: ZS ZEN 3600) at pH conditions repre- 
sentative of deep saline aquifers (6—12 pH) [62]. A greater ab- 
solute zeta-potential value (>30 mV) suggests that repulsive 
forces are dominant among nanoparticles, whereas a lower 
absolute zeta-potential value (i.e., near the isoelectric point, 
which is 0 mV) represents the dominancy of attractive forces 
among nanoparticles. These greater attractive forces result in 
a greater agglomeration and sedimentation ability of nano- 
particles. Generally, greater agglomeration is observed when 
particles are exposed to a saline environment [63—65]. 

The zeta potential and hydrodynamic particle size at 
different pH values are presented in Table 4. The data indicate 
that nano-alumina has zero isoelectric points in deionized 
water (when the pH values range from 8 to 10). The zeta po- 
tential decreased, whereas the average particle size increased 
with the increasing concentration of nanoparticles. Table 4 
further suggests that the alumina nanofluids are promising 
candidates for reversing the wettability of the mica-gas-brine 
system into more hydrophilic conditions due to the wide 
range of zeta-potential values (48 to —40 mV) [59]. Specifically, 
low nanoparticle agglomeration (158 and 206 nm) and higher 
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Fig. 1 — Total organic carbon (TOC) of the aged and unaged mica substrate. The blue line represents mica, and the gray line 
indicates organic aged mica. (For interpretation of the references to color in this figure legend, the reader is referred to the 
Web version of this article.) 


(b) 


Fig. 2 — Scanning electron micrographs of (a) a clean mica surface and (b) a mica surface aged in stearic acid/n-decane 
solution at 107° mol/L concentration. 


Table 3 — Mean elemental values of pure and aged mica substrates. 


Organic acids After aging Before aging 
Si (wt%) C (wt%) O(wt%) Al(wt%) K(wt%) Si(wt%) C(wt%) O(wt%) Al(wt%) K (wt%) 


Lignoceric acid WG di 50416 478418 190410 7.7410 225410 32+05 464415 192+11 88:0.7 
Stearic acid 20.8413 48416 481417 192410 7.2+09 223412 32+06 461409 19.7408 87+0.7 
Lauric acid 21.0413 44417 482420 189410 76+07 222+09 31407 468420 196+09 83+0.7 
Hexanoic acid AM Oar il4as BOs BOsesAfy iS eeils) bse ihil Ailijesily AQarO7/ E Ai Iss 0G /£)se 05 
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Zeta potential and hydrodynamic particle size of alumina nanofluid concentrations at different pH values. 


pH Zeta potential (mV) Average particle size (nm) 
0.05 wt% 0.1 wt% 0.25 wt% 0.75 wt% 0.05 wt% 0.1 wt% 0.25 wt% 0.75 wt% 
6 48.35 48.35 46.5 40.6 158 158 167 207.5 
8 17.62 16.85 11.545 10.305 551.6 540 494.4 525 
10 —6.13 —10.65 —14.35 —18 1056 642.5 512.5 540 
12 —28.4 —29.46 —38.8 —39.75 206 202 192.5 235.9 


stability (48.35 and —28.4 mV) is noticed at 6 and 12 pH fora 
0.05 wt% nanofluid concentration, whereas a higher agglom- 
eration (551.6 and 1056 nm) and lower stability (17.62 and 
—6.13 nm) is noticed at pH values of 8 and 10 for a 0.05 wt% 
nanofluid concentration. A similar trend was also observed for 
other nanofluid concentrations. Thus, the 
adsorption of the alumina nanofluid on the mica surface 
occurred through the definite interaction of —OH ions with the 
positive charges of the mica substrates [59]. 


irreversible 


Aging process of organic-aged mica in alumina nanofluid 

The mica substrates (aged in organic acid or n-decane solu- 
tions) were vertically immersed in different concentrations of 
nano-alumina formulations (0.05—0.75 wt%) for 5 h at 323 K 
and ambient pressure. The samples were submerged 
perpendicularly to minimize gravitational settling of nano- 
alumina at the ratio of 5:1; that is, for every 5 g of nanofluid 
dispersion, 1 g of mica was used. The nanoparticles were 
permanently adsorbed onto the mica surface due to the 
chemical reaction or interaction of the positively charged 
hydroxyl] group of mica substrate with the negatively charged 
hydroxyl] group of nanofluid dispersions. 

This phenomenon altered the wettability of the initially 
hydrophobic mica surfaces (originally aged in organic acid or 
n-decane) to an intermediate-wetting state when aged in 
alumina nanofluids. Generally, organic acid or n-decane mica 
surfaces became CO,-wet due to the esterification of the 
substrates with organic acids via a condensation reaction. 
However, the Al2O0; nanoparticles were irreversibly adsorbed 
on the aged mica to restore the original wettability state 
before the organic acid adsorption [50,55,59]. 


Measurements of the contact angle of mica-H-brine and 
mica—CO—brine systems 


Advancing (fa) and receding (0,) contact angles of the rock-gas- 
brine systems were measured at reservoir conditions (5, 10, 15, 
and 20 MPa and 308 K, 323 K, 343 K for pure mica and 15 and 
25 MPa and 323 K for organic- and nano-aged mica) using the 
tilted plate method (17° tilted angle) [50]. A detailed descrip- 
tion of the equipment and procedure for the contact angle 
measurements with their sources and types was described in 
our previous studies [26,35,59]. Initially, the mica substrates 
(pure, organic-aged, and nano-aged) were inserted into the 
sample holder (at 17° tilted arrangement), and the cell was 
firmly closed. Then, H3 or CO, was gradually injected to reach 
the desired temperature (308 K—343 K) and pressure 
(0.1-25 MPa). After the cell was filled with H2 or CO3 at the 
desired operating temperature and pressure, equilibrated 
brine droplets (5 + 0.75 uL) were introduced onto the tilted 


mica surface with a high-precision pointer. Subsequently, 6a 
and 0, were filmed at the trailing and leading edges of the 
droplet with a high-performance video camera [35]. The con- 
tact angles on the images were measured with image analysis 
software (Image J). For every experimental condition, mea- 
surements were replicated thrice, and the mean values of 6, 
and 6, were computed [66]. The standard deviation of the ex- 
periments at 0.1—-15 MPa was +3°, whereas that of experi- 
ments at 20—25 MPa was +5°. 


Computation of the solid-fluid interfacial tension and 
equilibrium contact angle 


Comprehensive procedures for calculating the solid-fluid IFT 
have been reported in the literature [27,43—47]. In addition, 
the relevant computation procedure is provided concisely in 
this section. 

As the equilibrium contact angle (0g) at thermodynamic 
equilibrium conditions is not accessible by experiment [67], 
the below equation calculates ðg: 


(1) 


NacoSO, + NRCOSO 
6: = arccos( AEE £) 


Na HNR 


where 0a and 6g represent experimental advancing and 
receding brine contact angles: 


(2) 


aa sin? 6, 3 mica sin? ôr 3 
A \2— 3cos0a + cos?20, id 2 — 3cosOp + COS? Oz 


Combining Young's equation (Eq. (3)) and Neumann's 
equations of state (Eqs. (4) and (5)) 


cosh; = “E = Ys (3) 
Yra 

Yse = Ys. + Yre — 2V YsLYLG [2 — b(YsL — ne)’| (4) 

Ys. = Ysc + Vic — 2V YscYrG [2 —B(Yse — tha) | (5) 


obtains Eq. (6), as follows: 


cosh; =1 af] B(Ys. ne)”| (6) 
YLG 


where ysc, Ys, and yr are the solid-gas, solid-liquid, and 
liquid-gas IFTs, respectively, and 6 is the adjustable constant 
[27]. 

Further, in Eq. (6), 0g and yzg are input parameters, and the 
yıc data at various pressures and temperatures were collected 
from Ref. [68], as presented in Figure S1 (supplementary file), 
and £ are the fitting parameters. Further, to mimic in situ gas 
geo-storage scenarios, the brine was not saturated with any 
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gas prior to yı measurements, and finally, ysg was deter- 
mined using Eq. (4). 


Results and discussion 


Effect of temperature and pressure on equilibrium contact 
angles and rock-fluid IFT of pure mica-gas-brine systems 


As a low-permeability clay mineral, mica can represent a 
proxy of caprock, and its wettability has a great influence on 
the structural sealing potential during gas geo-storage [26,35]. 
The geo-storage formation temperature and pressure change 
with depth [69]; thus, we investigated the effects of tempera- 
ture and pressure on the equilibrium brine contact angle (6g), 
rock-gas IFT (ysc), and rock-brine IFT (ys;) for pure mica-gas- 
brine systems, as illustrated in Figs. 3—5. 


Equilibrium contact angle of the pure mica-gas-brine system 

The effects of temperature and pressure on the equilibrium 
brine contact angle (0g) of the pure mica-gas-brine system 
were investigated and presented in Fig. 3. Clearly, 0g increased 
with increased pressure and decreased temperature. For 
example, in the UHS scenario (i.e., the gas is H3), 0g increased 
from 20° to 47° as the pressure increased from 5 to 20 MPa at 
323 K, and 6; increased from 19° to 41° when the temperature 
reduced from 343 K to 308 K at 10 MPa. In the CGS scenario (i.e., 
the gas is CO2), 0g rose from 32° to 79° with rising pressure 
from 5 to 20 MPa at 323 K, and bg increased from 47° to 62° once 
the temperature reduced from 343 K to 308 K at 10 MPa, 
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Fig. 4 — Effect of pressure and temperature on the pure 
mica-brine interfacial tension (IFT). 
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mica-gas interfacial tension (IFT). 


indicating that mica wettability shifted from intermediate- 
wet to weakly water-wet [8,70]. 

Fundamentally, pressure effects on the wettability (contact 
angle) are attributed to the intermolecular force and gas 
density [69,71,72], whereas the effect of the temperature on 
wettability (contact angle) is due to the intermolecular force 
and interfacial energy [73]. Generally, CO, and H; densities 
increased with increased pressure and decreased temperature 
[25]. Consequently, the intermolecular interaction between 
the gases and solid (rock surface) is higher at high pressures 
and lower at high temperatures [25,74]. 

According to the Young-Laplace equation, a larger contact 
angle leads to smaller capillary forces [76,77]. Thus, the gas 
structural trapping capacity could be reduced, resulting in a 
higher risk of gas leakage [78]. Thus, the formation tempera- 
ture and pressure (i.e., depth) effects must be considered 
regarding the structural storage capability and containment 
safety for UHS and CGS projects [78]. However, at the same 
thermo-physical conditions, 0g for Hə was always lower than 
g for CO3. For example, at 20 MPa and 343 K, the 6; values for 
H2 and CO, were 37° and 72°, respectively because H3 has a 
smaller gas density than CO,; thus, weaker intermolecular 
forces exist between H, and the mica surface compared to CO, 
and the mica surface [71,79]. 


Pure mica-brine interfacial tension 

Brine density is almost constant at various pressures [25,80], 
resulting in a weak dependency of the intermolecular forces 
between the mica and brine on pressure (Fig. 4). The mica-brine 
IFT (y5,) remained constant as the pressure increased from 5 to 
20 MPa. However, ys, is sensitive to temperature change 
because increasing temperatures reduce the liquid density, as 
the energetic brine molecules occupy more space upward. For 
example, for Hz, ys, reduced from 58 to 41 mN/m as the tem- 
perature rose from 308 K to 343 K. For CO, raising the tem- 
perature from 308 K to 343 K dropped ys, from 10 to 5 mN/m. 
Likewise, ys, for H? was much larger than that for CO2. 


Pure mica-gas interfacial tension 

The pure mica-gas IFT (ysc) decreased with increasing pres- 
sure for both gases, as depicted in Fig. 5. For instance, at 323 K, 
by raising the pressure from 5 to 20 MPa, ysg for H2 dropped by 
19 mN/m (from 114 to 95 mN/m), and ysg for CO, reduced by 
36 mN/m (from 48 to 12 mN/m). These results are qualitatively 
consistent with the pressure influence on basalt-gas and 
quartz-gas IFT values [27,46]. 

In contrast, the effect of temperature on ysg depends on 
the gas type. Specifically, at 10 MPa, ysc for H2 reduced by 
9 mN/m (111—102 mN/m), whereas ysg for CO, increased by 
20 mN/m (14—34 mN/m), as the temperature was elevated 
from 308 K to 343 K. 

The decreasing mica-gas IFT with increasing pressure can 
be ascribed to the adhesion of the CO, or H, onto mica sur- 
faces due to the increasing cohesive energy density of H3 or 
CO, resulting in favorable mica-gas interactions enhancing 
the mica-H, or mica-CO, wettability [43]. Fig. 5 further con- 
firms that ysg for Hə was larger than ysg for CO, at similar 
thermo-physical conditions because of fewer van der Waals 
interactions between the mica substrates and Hz molecules 
with lower density compared to CO, [43]. For instance, at 343 K 
and 15 MPa, the ys, values for H, and CO, were computed as 96 
and 22 mN/m, respectively. 

In principle, for a particular gas type, the rock-gas inter- 
molecular forces and gas density increase with increasing 
pressure [71,79]; thus, ysg reduces correspondingly [46]. For 
several gas types, at a similar temperature and pressure state, 
the density of CO is higher than that of H2. Consequently, Ysg 
follows the opposite order: ysg for Hz is larger than ysg for CO. 
This finding is qualitatively harmonious with the latest re- 
ported results on the influence of gas density on basalt-COz, 
shale-CO., and quartz-H2 IFT values [27,29,30]. In contrast, a 
lower ysg represents higher contact angles (and, thus, hy- 
drophobic wettability), which represents lower gas storage 
capacities (e.g., 0 > 90° may cause structural leakage) [49,54]. 


Effect of organic acids, alkyl chain length, and pressure on 
equilibrium contact angles and rock-fluid IFT of organic- 
aged mica-gas-brine systems 


As organic matter adsorption and contamination prevail in 
the subsurface (e.g., deep saline aquifers or depleted hydro- 
carbon reservoirs) [81—84], the influence of organic matter 
adsorption on the interfacial properties of caprock at reservoir 
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conditions must be explored. Therefore, in this section, the 
effects of the concentrations and types of organic acid (as a 
model organic matter) on the interfacial characteristics of 
mica were investigated (Figs. 6—8). 


Equilibrium contact angle of the organic-aged mica-gas-brine 
system 

The equilibrium brine contact angle (6g) of the organic-acid- 
adsorbed mica-gas-brine system decreased with the 
decreasing concentration and carbon number of organic acid 
at constant thermos-physical conditions (Fig. 6). For instance, 
at a constant pressure (25 MPa) and temperature (323 K) 
representative of CGS conditions, as the concentration of 
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stearic acid increased from 107° to 107? mol/L, 6g increased 
from 74° to 117°, indicating that the mica wettability shifted 
from intermediate-wet to weakly CO2-wet [70]. Similarly, at 
the same thermo-physical UHS conditions, 0g increased from 
60° to 95°, indicating that the mica wettability shifted from 
weakly water-wet to intermediate-wet [26]. With a constant 
pressure of 25 MPa, the organic acid concentration (10-°mol/ 
L ) and increasing carbon number (6—24, i.e., from hexanoic to 
lignoceric acid), 0g increased from 78° to 99° for mica—CO». 
—brine systems. Likewise, a similar trend was observed 
regarding the influence of organic acid concentration and 
carbon number on g under UHS conditions, where 6g 
increased from 65° to 84°, as depicted in Fig. 6. 
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Fig. 6 — Effect of pressure and various organic acid concentrations on the equilibrium contact angle of the mica-gas-brine 
system. Readers are referred to Refs. [26,85] for advancing and receding contact angle data for organic-aged mica-H2-brine 


and to Ref. [35] for organic-aged mica—CO,—brine. 
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Fig. 7 — Effect of pressure and various organic acid concentrations on the mica-brine interfacial tension (IFT). 


Increasing the concentration of organic acid and the car- 
bon number resulted in stronger adsorption of organic acid 
molecules on mica surfaces [86]. Such strong adsorption of 
organic acid on the substrate surface and a high organic acid 
molecular coverage can significantly modify the wetting 
condition of the surface toward a less hydrophilic state 
[49,51,53,54,87,88]. A strong repulsion of polar H2O molecules 
occurs from the hydrophobic surface, whereas the nonpolar 
H2 and CO; approach the hydrophobic surface more easily and 
are easily attached to the mica surface. Thus, the response of 
the observed øg to the organic concentration and carbon 
number in Fig. 6 can be reasonably explained. 


Organic-aged mica-brine interfacial tension 

The results reveal that the brine density does not affect the 
mica-brine IFT (ysı) with the increase in pressure (Fig. 7). 
However, when the concentration of the organic acid and 
carbon number increased, the adsorbed organic acid (ys ,) also 
increased, as depicted in Fig. 7. Specifically, under UHS con- 
ditions (at 15—25 MPa and 323 K), with the increasing lignoceric 
acid concentration from 107? to 107? mol/L, yg, rose from 51 to 
60 mN/m. Similarly, at the same thermo-physical CGS condi- 
tions, ys, rose from 6 to 18 mN/m, when the lignoceric acid 
concentration increased from 10-2 to 10-2 mol/L. Whereas, at 
15—25 MPa and 323 K, the constant organic acid concentration 
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Fig. 8 — Effect of pressure and various organic acid concentrations on the mica-gas interfacial tension (IFT). 


(10-3mol/L ) and increasing carbon number (6—24, i.e., from 
hexanoic to lignoceric acid), ys, increased from 53 to 58 mN/m 
for mica-H2-brine systems. Likewise, a similar trend was 
observed regarding the influence of the organic acid concen- 
tration and carbon number on ys, under the same CGS con- 
ditions, where ys; increased from 12 to 15 mN/m, as depicted in 
Fig. 7. These results prove that the extent of organic matter 
adsorption is a crucial factor influencing the rock-fluid IFT. 

In line with the mechanism of organic acid adsorption in- 
fluence on the equilibrium contact angle, stronger organic 
acid adsorption (induced by a greater concentration of organic 
acid and carbon number) causes a higher energy barrier for 
spreading water on aged mica substrate. Consequently, yg, 
values were higher with a higher carbon number and an 
increasing organic acid concentration because the surfaces 


become less water-wet when the IFT between the mica and 
brine is higher [43]. 


Organic-aged mica-gas interfacial tension 

The mica-gas IFT (ysç) decreased with an increasing concen- 
tration of organic acid and carbon number, as depicted in Fig. 8. 
For example, at 323 K and 25 MPa, under UHS conditions, Ysg 
dropped to 62 from 84 mN/m with the increasing concentra- 
tion of lauric acid from 10-° to 107? mol/L. Similarly, at the 
same thermo-physical CGS conditions, ys. dropped to 7 from 
15 mN/m with the increasing concentration of lauric acid from 
107? to 107? mol/L. Whereas, at 25 MPa and 323 K, the constant 
organic acid concentration (10-*mol/L ) and increasing carbon 
number (6—24, i.e., from hexanoic to lignoceric acid), ys¢ 
reduced from 10 to 3 mN/m for mica—CO,—brine systems. 
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Fig. 9 — Effect of pressure and various alumina nanofluid concentrations on the equilibrium contact angle of the 


mica—CO,—brine system. Readers are referred to Ref. [59] for advancing and receding contact angle data regarding alumina 
aged mica—CO,—brine. 
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Fig. 10 — Effect of pressure and various alumina nanofluid concentrations on the mica-brine interfacial tension (IFT). 


Likewise, a similar trend was observed regarding the influence 71 mN/m, as depicted in Fig. 8. This is because stronger organic 
of the organic acid concentration and carbon number on ysg acid adsorption causes a lower energy barrier for gas molecules 
under the same UHS conditions, where ysg dropped to 47 from to approach the mica surface. 
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Influence of the alumina nanofluid concentration and 
pressure on equilibrium contact angles and the rock-fluid 
IFT of the nano-aged mica—CO2—brine system 


Nanofluids are widely used for various subsurface applica- 
tions, including altering the wettability of a gas-wet system 
into a water-wet system [59,63,64,89,90] and modifying the 
interfacial properties of hydrophobic caprock [50,55,91]. 
Threshold concentrations of nanofluids are crucial, where 
wetting characteristics and interfacial properties of caprock 
can be improved to implement CGS projects successfully. 
Therefore, this section presents the investigation of the effect 
of various concentrations of alumina nanofluid on the 
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interfacial and wetting characteristics of the mica—CO,—brine 
system (Figs. 9—11). 


Equilibrium contact angle of the alumina aged mica—CO2—brine 
system 

Generally, 0g decreased initially with an increasing alumina 
nanofluid (Al,O3) concentration (Fig. 9), but 6g later increased 
with an increasing nanofluid concentration, which occurs 
when the optimum concentration threshold is passed [50,59]. 
At realistic CGS conditions, the mica caprock was first aged in 
various organic acids (at a concentration of 107? mol/L) to 
build an initial baseline. Then, the organic-aged mica surfaces 
were aged in different concentrations of alumina nanofluid 
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Fig. 11 — Effect of pressure and various alumina nanofluid concentrations on the mica-CO, interfacial tension (IFT). 
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(0.05, 0.1, 0.25, and 0.75 wt%). At 25 MPa and 323 K with 
lignoceric aged mica substrates, 0g decreased from 85° to 71° 
with an increasing Al,O3 concentration from 0.05 to 0.25 wt%, 
but as the Al,O3 nanofluid concentration further increased 
from 0.25 to 0.75 wt%, øg increased from 71° to 80°, as depicted 
in Fig. 9. 

These results suggest that 0.25 wt% of alumina nanofluid is 
the optimum nanoparticle concentration ensuring a signifi- 
cant reduction in the contact angle. The principal mechanism 
for the modification of the mica surface wettability (reduction 
in contact angle , 6;) when aged in alumina nanofluid is 
ascribed to the permanent adsorption of Al,O3 nanoparticles 
onto the surface of mica [50,59]. Precisely, Al,O3 nanoparticles 
adsorbed evenly on the surface of the aged mica substrate as 
the concentration increased from 0.05 to 0.25 wt%. The exis- 
tence of —OH charges on nano-alumina dispersions plays a 
prominent role in ensuring the success of the irreversible 
adsorption of Al,O3 nanoparticles onto the positively charged 
mica surface (with the —OH group presence) [59], attracting 
water molecules to the nano-alumina treated surface, 
whereas CO, is repelled from the mica surface. This phe- 
nomenon is similar to the reported alteration of the negatively 
charged quartz surface wettability by nano-silica (with nega- 
tive charges) and has been identified as a critical reason for the 
alteration of hydrophobic into hydrophilic surfaces when 
aged in nanofluids [27,50,55,92]. 

However, the agglomeration of nanoparticles, resulting in 
the increased surface roughness of mica substrates, may occur 
at certain nanoparticle concentrations (0.75 wt% or upward). 
Cassie-Baxter's model revealed that the attraction of water to 
surfaces reduces when the surface roughness of such water- 
wet surfaces increases [27,93,94]. Hence, ôg increased as the 
nanofluid concentration increased from 0.25 to 0.75 wt% 
because higher energy is required to adhere water molecules to 
the surface compared to gas molecules [27]. Moreover, at higher 
nanofluid concentrations, the agglomerated particles reside in 
the bulk solution instead of irreversibly adsorbing onto mica 
surfaces, hindering the migration of surface-active species. 


Alumina nanofluid aged mica-brine interfacial tension 
Consistent with the literature study on the effect of silica 
nanofluid for quartz-brine IFT [27], the mica-brine IFT (y5,) 
also decreased as the initial alumina nanofluid concentration 
increased and then increased with the further increase in 
alumina nanofluid concentration (Fig. 10). For example, for the 
lignoceric acid aged mica, at 15—25 MPa and 323 K, ys, reduced 
from 10 to 7 mN/m as the Al,O3 nanofluid concentration rose 
from 0.05 to 0.25 wt%, suggesting that the attraction and 
spread of water on the mica surface become favorable at the 
optimum nanofluid concentration. 

However, ys, increased from 7 to 9 mN/m with the further 
increase in the AlO; nanofluid concentration from 0.25 to 
0.75 wt%, suggesting that the CO. wettability of mica became 
more favorable (gas molecules are more attracted to the mica 
surface), whereas the surface became less water-wet. 


Alumina nanofluid aged mica-COĻ interfacial tension 

In comparison, the mica-gas IFT (ysg) increased with the in- 
crease in the initial alumina nanofluid (Al,O3) concentration, 
whereas ysg decreased again with the further increase in the 


Al,O3 nanofluid concentration (Fig. 11). For example, at 15 MPa 
and 323 K with the lauric acid aged mica surface, ysg increased 
from 22.4 to 24.1 mN/m with the rising concentration of Al203 
nanofluid (0.05—0.25 wt%), whereas ys, reduced from 24.1 to 
23.2 mN/m with a further increase in the Al,O3 nanofluid 
concentration from 0.25 to 0.75 wt%. The ysg datasets 
confirmed the results from the equilibrium contact angle and 
Ys, that an optimum nano-alumina concentration exists for 
altering mica substrates contaminated with organic acids 
inherent in geo-storage formations into less hydrophobic 
conditions. These nanofluid concentration thresholds should 
be determined to assess storage-rock and caprock sealing 
potentials successfully and avoid containment uncertainties 
during UHS and CGS projects. 


Conclusion 


The wetting characteristics and rock-interfacial properties 
(including r, ys,, and ysc) of caprock at reservoir conditions 
are essential parameters that govern fluid distribution during 
the injection and storage of CO, and H; in storage rock and 
caprock and the subsequent withdrawal of H, at the peak 
demand period [49,51,70,75,78,95—98]. However, due to the 
technical constraints of measuring these three parameters 
through laboratory experiments, the rock-fluid interfacial 
data, particularly related to the H3 sealing potential of cap- 
rocks, are rarely reported in the literature [27]. Consequently, 
we calculated these three parameters for CO, and H, geo- 
storage conditions. Moreover, the influences of the concen- 
tration and carbon number of organic acid and nanofluids 
were explored to ascertain the extent of wetting phenomenon 
modifications by organic contamination inherent in geo- 
storage formations and the tendency of Al,O3 nanofluids to 
restore the original wetting state. 
The following conclusions are drawn from this study: 


1. The øg increases with increasing pressure, increasing the 
concentration of organic acid and carbon number but de- 
creases with increasing temperature, suggesting that high 
pressure and the existence of organic contamination are 
unfavorable for CO, and H, storage. In contrast, warmer 
reservoirs could be ideal storage sites. 

2. The ysg reduces with increasing values of pressure, organic 
acid concentration, and carbon number, confirming that 
the adhesion and spread of water molecules to the rock 
surface are minimal, whereas the gas molecules wet the 
mica surfaces at high pressure and in the presence of 
organic acids with a high carbon number. 

3. Further, #¢, ys,, and ysg do not demonstrate a monotonous 
relationship with the Al,O3 nanofluid concentration as 0g 
and ys, decrease, whereas ysg increases with an increasing 
nano-alumina concentration from 0.05 to 0.25 wt%. 
Conversely, 0g and yg, increase, whereas ys decreases 
with an increasing nano-alumina concentration from 0.25 
to 0.75 wt%, implying that the optimum Al,0; nanofluid 
concentration (0.25 wt%) is essential for improving the 
sealing potential of the caprock. 

4. Higher contact angles and lower solid-gas IFT were obtained 
for mica—CO,—brine systems compared to mica-H>-brine 
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systems at all investigated conditions, indicating that the H3 
wettability of caprock is always less than the CO, wettability 
of caprock at similar thermo-physical conditions. 


This study provides critical insight into the structural UHS 
and CGS trapping capacities of caprock and the associated 
containment safety. Further, the work also guides subsurface 
modifications to optimize gas geo-storage conditions. This 
study presents the first documented rock-fluid IFT data for 
caprock representative minerals regarding UHS in downhole 
conditions. These data aid in the correct prediction of the 
structural sealing potential of caprock and suggest the con- 
ditions for modifications of the rock wetting phenomenon 
during UHS and CGS projects. 
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